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^ Abstract 

Bottom baryons decaying to a J/V' meson and a hyperon are reconstructed using 
1.0 fb"^ of data collected in 2011 with the LHCb detector. Significant J/ipA, 
—7- J/ip S~ and i?^ — t- J/i/j f2~ signals are observed and the corresponding masses 
^ are measured to be 

O Af(ylO) = 5619.53 ±0.13 (stat) ±0.45 (syst) MeV/c2, 

O M{S^) = 5795.8 ±0.9 (stat) ±0.4 (syst) MeV/c^, 

M{[2-) = 6046.0 ±2.2 (stat) ±0.5 (syst) MeV/c^, 

<^ while the differences with respect to the A9 mass are 

m 



M{S-) - M(ylO) = 176.2 ±0.9 (stat) ±0.1 (syst) MeV/c^, 
j>: ^(^^b~) - ^(^°) = 426.4 ± 2.2 (stat) ± 0.4 (syst) MeV/c^. 

^ These are the most precise mass measurements of the A^, and baryons to 

H date. Averaging the above A^ mass measurement with that published by LHCb 

using 35 pb-i of data collected in 2010 yields M(yl°) = 5619.44 ± 0.13 (stat) ± 
0.38 (syst) MeV/c2. 
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Hadrons are systems bound by the strong interaction, described at the fundamental 
level by quantum chromodynamics (QCD). While QCD is well understood at high energy 
in the perturbative regime, low-energy phenomena such as the binding of quarks and 
gluons within hadrons are more difficult to predict. Several models and techniques, such as 
constituent quark models or lattice QCD calculations, attempt to reproduce the spectrum 
of the measured hadron masses (for a review, see Ref. yij). While the masses of all expected 
ground-state mesons are now well measured, baryon data are still sparse. In particular, 
only six out of the sixteen 6-baryon ground states predicted by the quark model have been 
observed so far [2] . A complete and reliable experimental mass spectrum would allow for 
precision tests of a variety of QCD models [3] . 

The mass measurement of the heaviest observed b baryon, the state with bss 
valence quark content, is of particular interest. While both the DO and CDF collaborations 
have claimed the observation of the f2f^ — ?► J/ipQ^ decay, the reported mass values, 
6165 ± 10 (stat) ± 13 (syst) MeV/c^ from DO [Ij and 6054.4 ± 6.8 (stat) ± 0.9 (syst) MeV/c^ 
from CDF 0, differ by more than 6 standard deviations. On the other hand, there is good 
agreement between the mass measurements of the (bsd) baryon, which has also been 
observed by DO [6] and CDF [7j in the — )• J/ip E~ mode and, more recently, by CDF ^ 
in the Sj^ — )■ S'^vr" mode. These measurements average to 5791.1 ± 2.2MeV/c^ [2]. 

This letter presents mass measurements of the weakly decaying {bud), and f]^ 
baryons using the decay modes — )■ J/il^A, — )■ J/ipS~ and i?^ — )■ J/ipQ~ (charge- 
conjugated modes are implied throughout). The mass differences with respect to the A^ 
mass are also reported. This analysis uses data corresponding to an integrated luminosity 
of 1.0 fb~^ and collected in pp collisions at a centre-of-mass energy of a/s = 7 TeV with 
the LHCb detector in 2011. 

The LHCb detector ^9J is a single-arm forward spectrometer covering the 
pseudorapidity range 2 < 77 < 5, designed for the study of particles containing 6 or c 
quarks. The detector includes a high precision tracking system consisting of a silicon-strip 
vertex detector surrounding the pp interaction region, a large-area silicon-strip detector 
located upstream of a dipole magnet with a bending power of about 4 Tm, and three sta- 
tions of silicon-strip detectors and straw drift tubes placed downstream. The combined 
tracking system has a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 
0.6% at lOOGeV/c, and an impact parameter resolution of 20 |am for tracks with high 
transverse momentum. Charged hadrons are identified using two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
and a hadronic calorimeter. Muons are identified by a system composed of alternating 
layers of iron and multiwire proportional chambers. The trigger consists of a hard- 
ware stage, based on information from the calorimeter and muon systems, followed by a 
software stage which applies a full event reconstruction. 

Precision mass measurements require the momenta of the final state particles to be 
determined accurately. Therefore, an important feature of this analysis is the calibration 
of the tracker response. This accounts for imperfect knowledge of the magnetic field and 
tracker alignment [TT] . In order to reduce these dominant contributions to the systematic 
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uncertainty, a two-step momentum calibration procedure is applied. Firstly, inclusive 
J/ip — )■ fi'^fi~ decays are used to account for the changes in the relative momentum 
scale between different data taking periods. Secondly, the absolute scale is derived from 

—7- J/ipK^ decays, taking the known masses [2] as references. In this procedure, the 
use of a J/ip mass constraint allows the momentum scale to be determined as a function 
of the track kinematics. The resulting calibration is checked with a variety of fully 
reconstructed decays listed in Fig. [TJ For each mode the mass distribution is modelled 
taking into account the effect of QED radiative corrections, resolution and background, 
and the mean mass value is determined. Following the procedure described in Ref. pT] . 
the deviation of the measured mass from the expected value is converted into an estimate 
of an average momentum scale bias independent of time and track kinematics. The bias 
is referred to as a, which is defined such that the measured mass becomes equal to the 
expected value if all particle momenta are multiplied by 1 — a. By definition, one expects 
a = to a good precision for the — t- J/ipK^ calibration mode. Figure [l] shows the 
resulting values of a. The nominal mass measurements are performed with a = and 
the largest value of |a| amongst the considered modes, 0.3 x 10"'^, is conservatively taken 
as the systematic uncertainty on the calibrated momentum scale. This uncertainty is 
somewhat larger than the 0.2 x 10~^ achieved with the 2010 data pT], due to changes in 
the alignment of the tracking devices. 

In this analysis, the 6-baryon decays — )■ J/ipA, S"^ — )■ J/ipS~ and — )■ J/il^ f2~ , 
followed by J/ip — )■ fJ^'^f^~, A — )■ p7i~ , S~ — Att~ and f]~ — )■ AK~, are reconstructed. The 
topology of these decays is characterised by the long-lived particles in the decay chain. 
The lifetime of weakly-decaying b baryons is ~1.5 ps and their decay vertex is expected 
to be separated from the primary pp interaction vertex by ~6 mm on average. The J/ijj 
candidates are reconstructed from pairs of oppositely-charged tracks, originating from the 
secondary vertex, that have hits in the muon detector. In the i^b) decay chain the 
long-lived S~ decays into a A and a charged pion (kaon) at a tertiary vertex and 

the A decays at a quaternary vertex. Since ~90% of the decays are not fully contained in 
the vertex detector, tracks that have no hits in the vertex detector are also considered in 
the reconstruction of the tertiary and quaternary vertices. 

The selections of and candidates are identical apart from the choice of the 
S~, , i7~, f]^ invariant mass ranges and particle identification requirements on the 
pion (kaon) from the ) vertex. The A^ selection is slightly different, owing to the 

different topology. 

The J/ip candidates are required to satisfy |M^^ — Mj/^\ < A.2a where M^^ is the 
reconstructed di-muon mass, Mjj^ the J /ip mass [2] and a the estimated event-by-event 
uncertainty on M^^ (typically lOMeV/c^). The invariant mass windows for the A, S~ 
and i7~ candidates are ±6MeV/c^, ±llMeV/c^ and ±llMeV/c^ around the expected 
masses [2] , respectively. Particle identification requirements are applied to the kaon from 
the i7~ candidate and the proton from the A decay to improve the purity of the selected 
daughter particles, but none is placed on the pion from the S~ candidate. In addition, 
the hyperon decay vertices are required to be downstream of the 6-hadron decay vertex. 

The A^, and fl^ mass resolutions are improved by performing a fit of the decay 
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Figure 1: Average momentum scale bias a determined from the reconstructed mass of various 
decay modes after the momentum calibration procedure. The decays are divided into two 
categories according to whether both daughter tracks (a) have hits or (b) do not have hits 
in the vertex detector. The black error bars represent the statistical uncertainty whilst the 
(yellow) filled areas also include contributions to the systematic uncertainty from the fitting 
procedure, the effect of QED radiative corrections, and the uncertainty on the mass of the 
decaying meson [2]. The (red) dashed lines show the assigned uncertainty of ±0.3 x 10~^ on the 
momentum scale. 

topology and vertices [12] while constraining the masses of the J/ip , A, S~ and f2~ hadrons 
to have their known values [2], the final-state and intermediate long-lived particles to 
originate from common vertices according to the decay chain, and the b baryon to originate 
from the primary vertex. Three additional variables are considered for the selection. These 
are the P^r degree of freedom (x^/ndf) from the fit, the reconstructed decay time and 
the Xip of the b baryon from the primary vertex. The Xip is defined as the difference in the 
of the primary vertex fit with and without the 6-baryon candidate. In the case of the 
candidates, the selection requirements for these variables are chosen to maximise the 
expected significance of the signal; the same selection is used for the f2^ candidates. To 
determine the significance for a set of selection criteria the background yield is estimated 
from the yield of candidates found in mass side-bands in the ranges 5600-5700 MeV/ 
and 5900-6100 MeV/ c^. The expected signal yield is estimated using the product of the 
world average hadronisation fraction for b — )■ Sj^ and branching fractions for — )■ 
J/ipS' and subsequent daughter particle decays p], the bb production cross-section in 
the LHCb acceptance [13] and the selection efficiencies obtained from simulation. The 
selection criteria giving the highest expected signal significance correspond to a decay time 
greater than 0.25 ps, a x^/ndf smaller than 4 and a Xip smaller than 16. Amongst these, 
the decay time requirement is the most powerful given the high level of background close 
to the interaction point. In the case of the yl° candidates, the decay time is required to be 



3 



1800 (a) A[; ^ J/V|/ A 
1600 E- 

1400 
1200 
1000 

800 E- 

600 

400 

200 

§500 5550 



LHCb 



5600 5650 5700 5750 

M(J/\|fA) [MeV/c'] 




5900 6000 
M(J/vS") [MeV/c2] 



" 24 

:^ 22 

■Z 20 
% 18 

16 

Cii4 

» 12 

S 10 

•S 8 

1 ' 



(c) Q.^ J/\|/ £2 




LHCb 



5800 5900 6000 6100 6200 6300 

M(J/\|fJl") [MeV/c2] 

Figure 2: Invariant mass distribution for (a) — > J/i^A, (b) — > J/'ipS~ and (c) J?^ — > 
J/")/; 17~ candidates. The results of the unbinned maximum hkehhood fits are shown with sohd 
fines. 
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m 

maximum likelihood fit. The A^, , fl^ 



greater than 0.3 ps and the x^/iidf smaller than 5 (no requirement on the Xip is made). 
The possibility of a cross-feed background between and is investigated using 
simulation and found to be negligible in comparison with the combinatorial background. 
The invariant mass distributions of the selected A^, and i7j~ candidates are shown 
In each case, the mass is measured by performing an unbinned extended 

candidates are retained for the mass fit if their 
invariant mass lies in the range 5500-5750, 5600-6000, 5800-6300 MeV/c^, respectively. 
The signal component is described with a single Gaussian function (or the sum of two 
Gaussian functions with common mean in the case of the A'^ baryon) and the background 
is modelled with an exponential function. The widths of the A^ and signals are left 
unconstrained in the fit. Due to the low expected yield for the signal, the width of 
the Gaussian function describing the signal is fixed to the measured Sj^ signal width 
multiplied by the ratio of and widths from the simulation (8.2 MeV/c^ for and 
8.9 MeV/c^ for S^). The fit results are given m Table [ll 

The statistical significance of the f2^ signal is determined using simulated pseudo- 
experiments with background only. We determine the probability that, anywhere in the 
mass range between 5800 and 6300 MeV/c^, a peak appears with the expected width and 
a yield at least as large as that observed in the data. This probability corresponds to 6 
standard deviations, which we interpret as the statistical significance of the i?^ signal. 
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Table 1: Results of the fits to the invariant mass distributions. The quoted uncertainties are 
statistical. The signal is described by a double Gaussian function with widths ai and CT2; the 
fraction of the yield described by the first component is 0.58 ± 0.11. 





Signal yield 


Mass [MeV/c2] 


Width(s) [MeV/c2] 




6870 ± 110 


5619.53 ±0.13 


Gi = 6.4 ±0.5 
(T2 = 12.5 ± 1.3 


■^b 


111 ± 12 


5795.8 ±0.9 


7.8 ±0.7 




19 ±5 


6046.0 ±2.2 


7.2 (fixed) 



The systematic uncertainties are evaluated by repeating the complete analysis (in- 
cluding the track fit and the momentum scale calibration when needed), varying in turn 
within its uncertainty each parameter to which the mass determination is sensitive. The 
observed changes in the central values of the fitted masses relative to the nominal results 
are then assigned as systematic uncertainties and summed in quadrature. The systematic 
uncertainties are summarised in Table [2l 

The dominant systematic uncertainty is due to the momentum scale calibration de- 
scribed previously, which is assigned an uncertainty of ±0.3 x 10~^. A significant contribu- 
tion to this uncertainty comes from the overall detector length scale along the beam axis, 
which is known to a relative precision of 10~^ [H]. This translates into a ±0.13 x 10~^ un- 
certainty on the momentum scale, and is included in the overall ±0.3 x 10~^ uncertainty. 
Most of the uncertainty related to the momentum scale is removed in the measurements 
of the mass differences. 

The uncertainty on the amount of material assumed in the track reconstruction for 
the energy loss (dE/dx) correction has been found to be small [H]. It translates into an 
uncertainty on the /1° mass of 0.09MeV/c^, which we apply to all masses. 

The invariant mass of candidates is computed assuming the central value of the 
i7~ world-average mass [2]. The uncertainty of ±0.29MeV/c^ on this value is propagated 
as a systematic uncertainty. A similar, but smaller, uncertainty is estimated for the Sj^ 
and masses from the imperfect knowledge of the S~ and A masses, respectively. 

Two alternative fits for the A^ signal are performed: a first fit where the candidates are 
split into two categories depending on whether the daughter tracks have vertex detector 
information or not, each category being described with a single Gaussian function where 
the two Gaussian functions have a common mean, and a second fit using the sum of 
two Crystal Ball functions [15] with common peak value and otherwise unconstrained 
parameters. The second fit allows to take into account possible QED radiative corrections. 

The mass fit is repeated using as an alternative model either the sum of two 
Gaussian functions with a common mean, or a single Crystal Ball function. In the 
mass fit, the fixed Gaussian width is varied within both the uncertainty of the fitted Sj^ 
width and the statistical uncertainty of the width ratio from simulation. 

An alternative background model assuming a linear shape leads to negligible changes. 
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Table 2: Systematic uncertainties (in MeV/c^) on the mass measurements and their differences. 
The total systematic uncertainty is obtained from adding all uncertainties in quadrature. 



Source 






% 


^b 


^b-A'b 


Momentum scale 


0.43 


0.43 


0.31 


0.01 


0.12 


dE/dx correction 


0.09 


0.09 


0.09 


0.01 


0.01 


Hyperon mass 


0.01 


0.07 


0.25 


0.07 


0.25 


Signal model 


0.07 


0.01 


0.24 


0.07 


0.25 


Background model 


0.01 


0.01 


0.02 


0.01 


0.02 


Total 


0.45 


0.45 


0.47 


0.10 


0.37 



We also repeat the and mass fits in a restricted mass range of 5650-5950 MeV/c^ 
and 5900-6200 MeV/c^, respectively, and assign the resulting change as a systematic un- 
certainty. 

In summary, the and baryons are observed in the — )■ J/ipA, — 

J/iIjS~ and — )■ J/ipQ^ decay modes using 1.0 fb~^ of pp collisions collected in 2011 
at a centre-of-mass energy of ^/s = 7 TeV. The statistical significance of the observed 
f2f^ — 7- J/tp fl~ signal is 6 standard deviations. The masses of the h baryons are measured 
to be 

M(yl^) = 5619.53 ± 0.13 ± 0.45 MeV/c^ 
M{E^) = 5795.8 ±0.9 ±0.4 MeV/c^ 
M{f2^) = 6046.0 ±2.2 ±0.5 MeV/c^ 

where the first (second) quoted uncertainty is statistical (systematic). The dominant 
systematic uncertainty, due to the knowledge of the momentum scale, partially cancels in 
mass differences. We obtain 

M{S^)-M{A°^) = 176.2 ± 0.9 ± 0.1 MeV/c^ 
M{f2;) - M{Al) = 426.4 ± 2.2 ± 0.4 MeV/c^. 

A measurement of the A^ mass based on the 2010 data sample, M(yl|J) = 5619.19 ± 
0.70 ± 0.30 MeV/c^, has been previously reported by LHCb [H]. Since the new alignment 
and momentum calibration procedures differ from those applied in the previous study, a 
possible correlation between the systematic uncertainties related to the momentum scale 
can be neglected. Considering that the only correlated systematic uncertainties are those 
due to energy loss correction and mass fitting, the weighted average of the two yl^ mass 
measurements that minimizes the total uncertainty is 

M{Al) = 5619.44 ± 0.13 ± 0.38 MeV/c^ . 

These A^, and mass measurements are the most precise to date. They are 
compared in Table [3] with the single most precise measurements from ATLAS, CDF and 
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Table 3: Comparison of the 6-baryon mass measurements using the full 2011 data sample with 
the single most precise results from the ATLAS [16] , CDF [5l [T7j and DO [31, 16] collaborations, 
and with the PDG averages [2J. The PDG averages contain the results from CDF and DO as well 
as the measurement from LHCb performed with the 2010 data sample. The quoted errors 
include statistical and systematic uncertainties. All values are in MeV/c^. 









Mm 


ATLAS 


5619.7 ± 1.3 






CDF 


5619.7 ± 1.7 


5790.9 ±2.7 


6054.4 ±6.9 


DO 




5774 ± 19 


6165 ± 16 


PDG 


5619.4 ±0.7 


5791.1 ±2.2 


6071 ± 40 


LHCb 


5619.5 ±0.5 


5795.8 ± 1.0 


6046.0 ± 2.3 



DO, and with the current world averages [2]. The A^ and results are in agreement with 
previous measurements. The result is in agreement with the CDF measurement [S], 
but in disagreement with the DO measurement 
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